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Abstract 



We present a new generating algorithm to construct exact non static so- 
lutions of the Einstein field equations with two-dimensional inhomogeneity. 
Infinite dimensional families of G\ inhomogeneous solutions with a self in- 
teracting scalar field, or alternatively with perfect fluid, can be constructed 
using this algorithm. Some families of solutions and the applications of the 
algorithm are discussed. 



Recently there has been a considerable interest in scalar field cosmologies. 
This interest is related to the studies of inflationary models where the scalar 
field (inflaton) acts as the main driving force for inflation. On the other hand, 
one of the most important issues related to the inflationary scenario is how 
sensitive are the inflationary models to the change in the initial conditions 
[[[]]. The dependence of the inflationary scenario on the initial conditions is 
usually studied either numerically or perturbatively, under some simplifying 
assumptions, or analytically, within the framework of exact solutions using 
cosmological models with a high degree of symmetry. Until now, it was only 
possible to study exact scalar field cosmological models with one dimensional 
inhomogeneity [@]. 

The main purpose of this letter is to make a step further in the study of 
exact inhomogeneous cosmologies by presenting a new generating algorithm 
to obtain solutions with a single isometry. The algorithm allows one to 
obtain spacetimes with homogeneity broken in two directions from already 
known ones with one dimensional inhomogeneity. One starts with a vacuum 
orthogonal G2 cosmology which is well studied and some general solutions 
for certain topologies are known (for a recent review see Ref. ||.) Then, by 
using a known algorithm the solution is transformed into a massless scalar 
field G2 cosmology. The second step involves a conformal transformation 
which breaks the G2 symmetry leaving one with the solution with the G\ 
symmetry but with a self interacting scalar field with a Liouville type of 
potential. 

Our algorithm generalizes a recently proposed solution generating tech- 
nique developed by Fonarev |4j], which is only valid for static metrics. 

To this end we consider a vacuum G2 cosmology with the following line 
element 

ds 2 = e /(M) (-dt 2 + dz 2 ) + G(t, z) (e p{t ' z) dx 2 + e - p{t ' z) dy 2 ) . (1) 

If p v , G v and f v represent a vacuum solution of the Einstein field equations 
then, by using a simplified version of the Charach-Malin || algorithm the 
following solution 

p = Bp v + C\ogG v 

f = f v + E Pv + F\ogG v (2) 
4> = A Pv + D\ogG v 
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will be a solution of Einstein-masless scalar field equations, where the scalar 
field 4> is minimally coupled to gravity and its energy momentum tensor is 
given by 

T a p = <j>, a $,p - - g a /3 0, 7 0' 7 • (3) 
The constants A, B, C, D, E and F are subject to the following constraints 

BC - E+2AD = 

C 2 + 2D 2 -2F = (4) 
B 2 + 2A 2 = 1 

Theorem: Given a minimally coupled massless scalar field solution g a } } de- 
scribed by the equations (l)-(3), a new solution with minimally coupled 
scalar field with exponential potential can be built by the following confor- 
mal transformation 

g ab = tt(x) 2 g ab . (5) 
The new scalar field (f> will be given by 

= k log ft (a;) (6) 

and the potential will take the form 

V(<f>) = V e-^. (7) 

For k 2 = 2 the conformal factor Q is given by ft = exp(aa;) and Vq = 2a. 
For any other value of k the conformal factor is given by 

2 12 — 2k 2 

A = V Q = {k2 _ 2)2 a 2 (8) 

The constants A, B, C, F and E can be expressed in terms of the constants 
D and k which are related with the dynamical and potential parts of the 
scalar field respectively: 



A = ± 
B = ± 



VWT2 
k 



VWT2 
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(k(Dk - 1) + 2D) 



(9) 



We will not prove the theorem here, but the results can be verified in 
straightforward manner by direct substitution into the Einstein field equa- 
tions. 

Before going any further several remarks are in order. This theorem gives 
one a simple way to reduce the symmetry of the problem. Generically, one 
reduces the symmetry at the cost of introducing a potential term. However 
in the special case k 2 = 6 in the expression (8) one still obtains solutions 
with a massless scalar field (or stiff fluid.) One does not have to start with 
the G2 symmetry but with any larger group which contains the G 2 - For 
example, one can start with G3 Bianchi I model and obtain from it the G2 
inhomogeneous solutions. Also, since we have not specified the character of 
the transitivity surface area which gradient GpG^ can be timelike, space- 
like or null one can construct models starting from the metrics describing 
cylindrical gravitational waves, plane waves, cosmic strings, etc.. The new 
obtained metrics can be looked at as perfect fluid spacetimes, if one wishes, 
by identifying the scalar field with the velocity potential of the fluid. Thus 
this is the first algorithm to obtain perfect fluid solutions with G\ symmetry. 

To demonstrate the working power of the algorithm we present several 
examples. 

Bianchi I as a seed. We start with the G3 Bianchi I vacuum (Kasner) solution 
by specifying the line element (1) with the metric functions 



The massless scalar field solution will be described by the metric functions 



G v = t, p v = (3 logt, /, 



(3 2 - 1 



log* 



(10) 



V 
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(11) 



while the scalar field will be given by 



<P= (Afi + D) log*, 



(12) 
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where the constants are subject to conditions (4). And finally, the solu- 
tion with a self interacting scalar field with an exponential potential will be 
described by the expressions (5)-(7), where the conformal factor is 

2 

Q(x) = when k ^2 

Q(x) = e ax when k 2 = 2 (13) 

The new solution represents an inhomogeneous G2 cosmology. This solu- 
tion is a particular case of more general solutions obtained in integrating 
directly the Einstein equations. 

Cosmological gravitational waves as a seed. Here one starts with the vac- 
uum cosmological model describing propagation of inhomogeneities in form 
of gravitational radiation on a spatially flat homogeneous background. To 
simplify the discussion we concentrate on a single mode perturbation. One 
could easily start from either an infinite-dimensional Gowdy like vacuum cos- 
mology or inhomogeneous cosmologies filled with monochromatic waves and 
pulses This, however, will not contribute qualitatively to our discussion. 
The vacuum solution is given by 

G v = t 

p v = (3 logt + A cos(uz) Jo{ujt) 



2 -1 1 r 

fv = — - — log t + (3 A cos(ujz) J (cut) + -(cot) 2 Al [Jl(ujt) + J 2 (ut) 

-^utAl cos 2 (u}z) J {uit) Ji(ut) (14) 

From this solution one constructs the minimally massless scalar field so- 
lution and then using the theorem the exponential potential solution. Since 
the construction is straightforward we will not present the explicit form of 
the solution. The resulting new solution represents a spacetime with two di- 
mensional inhomogeneity. The interpretation of the solution as representing 
propagation of gravitational waves is not at all clear now. It would be inter- 
esting to consider if the obtained solution represents the interaction of the 
gravitational waves with the scalar field. Work is in progress in this direction 
and we hope to be able to present results in the future. 

The structure of the solution is much more complicated than that ob- 
tained from the Bianchi I seed. Since the new solutions are always confor- 
mally related with the seeds the asymptotic behaviour at t = or at late 
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times can be easily deduced from the asymptotic behaviour of the seed met- 
ric. For example, since usually the G2 cosmologies with massless scalar field 
at late times tend to the so called Doroshkevich-Zeldovich-Novikov spatially 
anisotropic universe || the new solutions with self interacting scalar field will 
represent an inhomogeneous generalization of DZN universes at late times. 
Near the initial singularity, on the other hand, the structure is more com- 
plicated. It is known that the conformal transformations usually change the 
nature of the singularity. 

Also, the presence of the potential term for the scalar field changes the 
energy balance which can result in violating the strong energy condition. It is 
known that the breaking of the strong energy condition is intimately related 
with inflation. Note, that in a highly inhomogeneous model it is difficult, 
and sometimes even impossible to see whether the model inflates. If the 
cosmological model does not deviate strongly from homogeneity, say a slightly 
perturbed model, then the hypersurfaces of homogeneity perfectly define a 
preferred observer which can decide as to whether the model inflates. If the 
scalar field is not very inhomogeneous the hypersurfaces of cj) = constant 
serve to select again a preferred observer related to the motion of the perfect 
fluid which is defined from the scalar field. Thus, the notion of inflation in 
regular spacetimes is well clear. If the spacetime is highly irregular, on the 
other hand, neither the scalar field nor other geometric considerations help 
to indicate inflation. In this case one has to use somewhat a weaker way to 
specify inflationary behaviour. For example, one may look at the fulfillment 
of the strong energy condition the breaking of which is a necessary condition 
for a model to inflate. Obviously the models which satisfy the strong energy 
condition can be considered as the standard ones. We have checked the 
fulfillment of the strong energy condition for both models above. 

In the inhomogeneous models obtained from the Bianchi I seed (the mod- 
els with one dimensional inhomogeneity) the behaviour to respect inflation 
is quite simple. In general, at early times and for k 2 < 2 the strong energy 
condition is fulfilled and the cosmological model behaves in the standard way. 
Later the energy condition breaks down and the model starts to inflate. This 
also depends on the spatial coordinate x. In the models obtained from the 
G2 seed the behaviour presents much more complicated structure. In Fig. 
1 we have depicted the function 3p + p in terms of coordinates x and z for 
a constant value of time t. One can easily see that there is a "wave-like" 
pattern along the z direction, with energy condition broken periodically. As 
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the time passes the model starts to inflate everywhere, and the "wave-like" 
pattern gets washed away. 

Two solitons as a seed. One may start with a vacuum cosmological model 
describing the propagation and interaction of two strong gravitational pulses 
in a homogeneous background 0. The solution for the functions p v &ind G v 
is given by: 

G v = t 

p v = piogt — d\ cosh" 1 — ^ — - — d 2 cosh" 1 — 2 (15) 

where the constants z\ and z 2 can be complex in which case the real part 
of the function p v must be taken. The vacuum seed behaves as follows |J: 
at early times the model is highly irregular while at late times the solution 
tends to the background which is of Kasner type. Applying the generating 
algorithm one obtains a model which will have two dimensional inhomogene- 
ity at early times while at late times the model will tend to those obtained 
from Bianchi I seed. 

In the case when the soliton "poles" are real {z\ and z 2 are real) it is 
well known that the former seed solution represents the interaction region 
of the collision of two plane gravitational waves on a flat background. It is 
tempting, therefore, to interpret the new solution obtained after applying 
the algorithm as a collision of gravitational waves in the presence of a self 
interacting scalar field. This solution can give one some ideas as to how 
gravitational waves interact in the presence of a non trivial fluids. 

To this end it is important to stress that since there are a great amount 
of known G 2 solutions with different physical interpretations, topologies, sin- 
gularity structures etc. @ the new algorithm represents a powerful tool to 
obtain new families of G\ solutions with two dimensional inhomogeneities. In 
certain cases, imposing additional conditions on the scalar field (the timelike 
character of its gradient) the new solutions can be considered as perfect fluid 
ones. 

We have tried to present in this Letter a brief insight on the possibilities 
of the application of the new algorithm. We believe that the algorithm 
presented here opens a new door into the study of exact solutions of Einstein 
equations depending on three variables. The exact solutions obtained with 
this algorithm could also serve as test solutions for numerical relativity. The 
analysis and the study of various particular cases are left for future works. 
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Figure 1: The behaviour of 3p + p in terms of coordinates x and z at t — 2, 
when k 2 = 1, /5 = and D = o; = Ao = a; = l. The range of the coordinates 
x and z is -4 < x < 4 and -10 < z < 10. 
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